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Sample Calculations of Gamma-Ray Penetration 

Into Shelters: 

Contributions of Sky Shine and Roof Contamination 

Martin J. Berger and James C. Lamkin 

An approximate method is presented for calculating the penetration of gamma radia- 
tion in shelters. Sample calculations, for an assumed source energy of 1 million electron 
volts, are given for the following problems: (1) Dose rate inside houses and underground 
shelters whose roofs are covered with radioactive fallout, and (2) dose rate in open holes 
due to reflected radiation (sky shine) from fallout contamination on the surrounding ground. 
A detailed examination is made of the dependence of the dose rate in a shelter on the shape 
of the shelter, and on the position of the detector within the shelter. The estimated accuracy 
of the calculations is ±30 percent. 



1. Introduction 

The shielding against gamma radiation provided 
by shelters, houses, and similar structures is of con- 
siderable practical interest. The experimental data 
on this problem are usually obtained under complex 
circumstances, so that they are difficult to interpret. 
It is desirable, therefore, to obtain some theoretical 
understanding of the shielding problem. Calcula- 
tions with the required accuracy (say 30 to 50%) 
appear feasible, and are perhaps easier than experi- 
mental determinations. 

The well-developed semianalytic theory of the 
penetration and diffusion of gamma radiation 
[1,3] 2 treats only infinite homogeneous media, and 
is thus not directly applicable to the problem under 
consideration which involves more complicated 
boundary conditions. Yet the theory provides the 
foundations on which the calculations can be carried 
out. It would be possible, but rather laborious, 
to take the boundary conditions into account by 
means of random sampling (Monte Carlo). An- 
other approach appears more promising — the adapta- 
tion of infinite-medium results to realistic conditions 
through a suitable schematization. 

It is the purpose of this paper to explore one such 
schematization that can be applied to calculations 
of shielding against fallout radiation. In this schemat- 
ization, the radiation dose inside a shelter is com- 
puted as a suitable weighted integral over the angular 
distribution of the radiation dose in an infinite 
medium. The present paper contains only pilot 
calculations. More extensive and systematic cal- 
culations are now in progress as part of a National 
Bureau of Standards program of shelter evaluation, 
carried out under the auspices of the Federal Civil 
Defense Administration. 



i Work supported by the Federal Civil Defense Administration. 

2 Figures in brackets indicate the literature references at the end of this paper. 



2. Schematization 

The nature of the schematization can be most 
easily explained by applications to specific problems. 
Generalization to other situations is then straight- 
forward. Figure 1 illustrates the schematization 
to be employed for three typical situations that will 
form the basis for our discussion. 

2.1. Underground Shelter With Roof Covered by 
Fallout 

The shelter roof and the surrounding ground are 
assumed to be covered with radioactive material 
that constitutes a plane isotropic source of gamma 
radiation of constant strength per unit area, which 
emits gamma rays uniformity in all directions. 

Simplification 1 . Different materials of low aver- 
age atomic number, such as air, concrete, earth, or 
water, have nearly the same gamma-ray mass 
attenuation coefficients. Provided that all distances 
from the source are expressed in units of mass per 
unit area, it is therefore a good approximation to 
replace the combination air-concrete (or earth) by a 
concrete medium that is infinite and homogeneous, 
except that it contains an air cavity (the inside of 
the shelter). 

Simplification 2. In order to reach a detector 
somewhere inside the shelter, radiation must pene- 
trate varying amounts of concrete or earth, depend- 
ing on its path. Radiation coming through the 
ceiling has the least amount of material to traverse 
and is therefore least attenuated. There are two 
types of paths on which photons could reach the 
detector through the floor. They can detour the 
shelter to come up from below, but in this case the 
attenuation is so strong that the contribution of the 
detoured photons is negligible. There is also the more 
important possibility for photons to enter the shel- 
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Figure 1. Schematization. 



ter through the ceiling, to penetrate into the 
floor (or side walls) and then to be backscattered. 
Some radiation will also penetrate to the detector 
through the side walls of the shelter. But because 
of the small size of a mean free patli in concrete or 
earth (a few inches) only a small side-wall region 
near the ceiling can make a nonnegligible contribu- 
tion, resulting in an effective ceiling area somewhat 
larger than the true ceiling area. Both the back- 
scattering correction and the side-wall penetration 
correction are small. They are briefly discussed in 
sections 6.1 and 6.2 but are omitted from the 
initial schematization. It will be a good first ap- 
proximation to consider only radiation entering 
through the ceiling*. The mean free path in air is 
of the order of several hundred feet, and is very 



large compared to the shelter dimensions. The 
radiation thus travels in a straight line inside the 
shelter and can reach a detector only if its direction 
of motion lies within the solid angle subtended by 
the ceiling with respect to the detector. This solid 
angle is indicated in figure 1 by the label "allowed 
directions. " The notion of "allowed directions" will 
be formalized later through the introduction of an 
"angular response function." 

Summary of the Calculation. The calculation con- 
sists of the following steps: (1) Replacement of the 
air-shelter complex by an infinite concrete medium; 
(2) computation of the amount and angular distribu- 
tion of the radiation in such a medium, at a distance 
from an infinite plane isotropic source equal to the 
thickness of the shelter roof; (3) computation of the 
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angular response function for a specific shelter- 
detector configuration; (4) computation of the total 
amount of radiation received by the detector, by an 
integral over the angular distribution weighted by 
the angular response function; and (5) corrections to 
account for additional radiation reaching the de- 
tector through the side walls, and by the reflection 
from interior surfaces of the shelter. 

2.2. House With Roof Covered by Fallout 

As indicated in figure 1, the schematization for 
this problem is essentially the same as for the 
previous problem. A number of points require com- 
ment, however. 

(1) The contamination now is confined to the roof 
so that the plane-isotropic source is Unite. Yet the 
situation is practically equivalent to that of an 
underground shelter underneath an infinite con- 
taminated plane. Note that in the underground 
shelter, fallout from regions other than that directly 
above the roof does not contribute appreciably to 
the dose, because of the large amount of intervening 
earth. 

The equivalence is, of course, only approximate. 
To illustrate the order of magnitude of the error in- 
volved, one may consider the following point. In 
the schematization, approximately 10 to 15 percent 
of the energy emitted by the source is reflected back 
toward the roof. In the actual situation, with air 
above the roof, this reflected radiation would spread 
laterally, and most of it would come down on the 
ground surrounding the roof. The reflected radiation 
assumed incorrectly to be incident on the roof is low 
in energy and diffuse in direction, so that its pene- 
trating power is much less than that of the source 
radiation. Hence, the resulting dose-overestimate 
will only be a few percent. 

(2) There is a somewhat better chance than in the 
case of an underground shelter, that radiation origi- 
nating on the roof will reach the detector through the 
side walls. But this radiation consists mainly of 
photons which after repeated air scattering are inci- 
dent on the side walls with low energies and therefore 
are not very penetrating. 

(3) If the ground surrounding the house is also 
covered with fallout, the penetration of radiation 
through the side walls becomes very important. The 
numerical details of this case are not worked out in 
the present paper, but the same schematization could 
be applied as in the case of a roof contamination, 
through the introduction of appropriate allowed 
directions and angular response functions defined 
with respect to the side walls, or perhaps the windows. 

2.3. Open Hole Surrounded by Fallout 

Again, the schematization is essentially the same 
as in 2.1. The detector is now in a position where it 
cannot "see" the source, i. e., only scattered radiation 
can reach it. This scattered radiation will consist 
mainly of photons that have gone up into the air and 
are reflected into the hole by scattering. This radia- 
tion component is customarily called sky shine. The 



radiation (list ributions that must now be used in steps 
(2) and (4) of the schematized calculation are those 
pertaining to the scattered flux at the source plane. 

3. Mathematical Formulation 

An isotropic source of gramma radiation is as- 
sumed to be located in the plane 2=0. Let I D (z, 6) 
sin Odd denote the dose delivered to a detector in an 
infinite homogeneous medium at a distance z from 
the source plane, by photons incident on the detector 
at angles between 6 and 6-\-dd (with respect to the 
z-axis) . Dose means the quantity technically called 
"absorbed air dose", i. e., energy dissipated in air, 
measured in rads (100 ergs/g) [2]. It will be con- 
venient to express the dose distribution in the form 



i D {z,e)=f{z)g{z,e), 



(i) 



where g(z,d) is an angular distribution normalized to 
unity : 



J g(z,6) sin 6d$=l. 



(2) 



The factor j(z) will represent the attenuation of 
radiation by the roof of thickness z in the problems 
of the house or underground shelter. For the open 
hole, one must replace j(z) and g(z,0) by quantities 
j s (z) and g s (z,B) pertaining to scattered radiation. 

Let the shelter-detector configuration be such that 
only a fraction \f/(6) of the photons with direction 
can reach the detector. This function is called the 
angular response function and is further discussed 
in section 5. The dose D(z) received by the partially 
shielded detector in the shelter is expressible as the 
product of three factors: 



where 



D(z)=Cf{z) G(z), 
G(z)= ' g(z,d)ip{0) am Odd, 



(3) 



(4) 



and C is a correction factor taking into account the 
contribution of radiation penetrating through the 
side walls or reflected from the interior surfaces. 
The determination of the geometrical protection 
factor G is the main objective of this paper. 

4. Calculation of the Radiation Field 

The standard solution for an infinite homogeneous 
medium is provided by the moment method [1]. 
Tabulations based on this method are available 
from which one can extract information about the 
dose f(z) from a plane isotropic source [3}. 3 The 
calculation of the angular distribution g(zfl) has been 
described in [4]. 

The numerical computations of this paper all per- 
tain to a source emitting 1-Mev photons. This 
source energy is a little higher than the average 
energy of fallout radiation, but close to the energy 

3 As part of the work sponsored by FCDA, further tabulations, with applica- 
tions to shelter evaluation, are now being prepared at NBS. 
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Figure 2. Angular distribution of absorbed air dose from a 
1-Mev plane isotropic source in an infinite medium. 

The calculation is for a water medium, but is in very good approximation also 
applicable to concrete. 

(a) g{z,d), at a distance noz=4 from the source plane. 

(c) g a (z,6), at the source plane wz=Q (This distribution represents only the 
contribution of scattered photons.) 

The curve (b) represents an isotropic distribution. 



of Co 60 -radiation, which is most suitable for con- 
trolled experimentation. The choice was moti- 
vated by the ready availability of results for g{zfi). 
In any case, calculations with a 1-Mev source will 
be sufficient for the illustrative purposes of this 
paper. 

Let j(z) be the dose in an infinite concrete medium 
at a distance z from a plane isotropic source. The 
results of a calculation by the moment method can 
be summarized by the approximate formula 



with 



/(*)=/o(s) +/.(*) 

fo(z)=kE 1 (ii z), 



(5) 
(5a) 



representing the contribution of unscattered radia- 
tion, and with 



f s (z) = l.04ke-°' 915 ^ z , 



(5b) 



representing the contribution of scattered radiation. 
The range of validity of (5b) is 0<mo^<8, E x is the 

exponential integral [5], and /x = 0.0635 cm 2 /# is the 
mass attenuation coefficient for gamma radiation in 
concrete at the source energy (1 Mev) . The constant 
k determines the source normalization. If the source 
plane is covered with 1 curie of gamma emitter per 
square foot, and the dose rate is expressed in rads 
per hour, the normalization constant has the value 



fc=32.2. 



(5c) 



If the fallout is spread over a rough surface, the 
schematization of the source as a perfect geometric 
plane is not altogether suitable. A possible alterna- 
tive schematization consists in assuming that the 
radioactive material is mixed homogeneously with a 
thin top layer of the ground of thickness t. In this 
case, the dose function f(z) must be replaced by 



f*V 



.«-y 



■z+t 



f(z')dz'. 



(6) 



Three angular distributions are shown in figure 2: 
(a) g(z, 6) at a distance n z=4 (0.897 lb of concrete 
per square inch) ; (b) a completely isotropic distri- 
bution; and (c) the distribution g s (z,6) for scattered 
radiation only, at the source plane fjL z—0. The 
distribution (a) peaks at 0=0° and is for use in the 
problems treated in sections 2.1 and 2.2. The dis- 
tribution (c) peaks at 0=90° and is for use in the 
problem in section 2.3. Both of these distributions 
pertain to an infinite plane-isotropic source and 
were calculated according to the method of reference 
[4], The distribution (b) is flat (when plotted against 
cos 6) and has been included as an intermediate case. 
Together, the three cases span a wide range of 
possible distributions. 

5. Calculation of the Response Function 

To each angle 6 there corresponds a family of unit 
vectors (sin0cos<p, sin0sin<£>, cos0), indicating possible 
photon directions. If these unit vectors are required 
by the schematization to lie within a specified solid 
angle, then not all values of the azimuth <p are 
permitted, but only azimuths in a region Ae<2ir. The 
response function is defined as the fraction ^(0) =Ao/2t. 

5.1. Cylindrical Geometry 

As a simple first example, we consider a shelter in 
the form of an upright cylinder of height H and 
radius R, with the detector located at the center of 
the bottom of the shelter (see fig. 3). The response 




Figure 3. 



Calculation of the response function for a cylindrical 
shelter. 
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function is in this case particularly simple: 



<6<d 1 =cos- 1 



H 



=,A*=2ir,f(0) = l 



0i<0<tt ,A e =0, ^(0)=O. 

5.2. Parallelepiped Geometry 



(7) 



We consider a box-shaped shelter of length L, 
width W, and height H. It can be assumed without 
loss of generality that L > W. The response func- 
tion depends on solid angles, and is thus a function 
of the shape but not of the absolute size of the 
shelter. Therefore, we eliminate one parameter by 
dividing all dimensions of the shelter by L, so that 
the scaled length is unity, the scaled width is 



W/L = e<l, 

and the scaled height is 

H/L=h. 



(8) 



(9) 



The construction of the response function is 
illustrated in figure 4 for a detector in a bottom corner 
of the shelter. 
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Figure 4. Calculation of the response function for a box- 
shaped shelter. 



Angular range 


Response function, \p(6,€,h) 


O<0<0 (1 — ros -1 


1 
1 • -i / e ^ 




V/* 2 +i 


2tt 3m \h tan 6 ) 


V^ 2 +€ 2 +l 

e 2 <e<TT 


2tt Dm \h tan e) 2ir C °" \h tan dj 




(10) 



The response function (10), evaluated with the 
parameters €=0.5 and h=0.5, is shown in figure 5. 

5.3. Combination Rule 

Expression (10) can readily be generalized so as 
to provide the response function for an arbitrary 
detector location. In the first place, it is to be noted 
that — within our schematization — the vertical dis- 
tance of the detector from the shelter ceiling (or top 
of the open hole) must be assumed as the effective 
height II. Now consider a detector located at the 
bottom of the shelter, in a horizontal position P with 
coordinates (a,0e), as illustrated in figure 6. We 
now imagine that the shelter is partitioned by two 
vertical walls (perpendicular to each other), which 
pass through the point P, and are very thin, but 
perfect, radiation shields. The detector can then 
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Figure 5. Response function for a detector at a bottom corner 
of a box-shaped shelter with scaled width (horizontal eccen- 
tricity) e=0.5 and height h = 0.5. 
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Figure 6. Diagram illustrating the application of the com- 
bination rule by which response functions for off-corner 
positions are calculated. 



be thought of as being located at the corner of each 
of four adjoining box-shaped shelters with the 
dimensions listed below. 



3 


e/ 


*f 


1 


min (<x, j8e) 
max (a, /3e) 

min (1 — a, /3e) 
max (1 — a, /3e) 

min (a, e— (3e) 
max (a, e— /3e) 

min (1 — a, e — jSe) 
max (1 — a, € — jSe) 


A 


2 

3 

4 


max (a, /3e) 


max (1 — a, /3e) 
ft 


max («, e — |8e) 
ft 


max (1 — a, e — jSe) 


0<«<1, 0</3<l. 

min (x, i/)= smaller of x and y. 

max (x t y) = larger of x and y. 



The response function \f/ P at point P is 



5.4. Legendre Polynomial Expansion 



(ii) 



In order to obtain the geometric protection factor 
G(z), it is necessary to evaluate the integral in (4) for 
each distribution g(z,d) of interest. When many dis- 
tributions must be included, this procedure tends to 
become laborious. It is then convenient to perform 
a set of standard integrations based on an expansion 
of the response function into a series of Legendre 
polynomials: 



i=o * 



(12) 



The theory of radiation diffusion provides as the 
most direct result, a similar Legendre potynomial 
expansion for the angular distribution: 

g(z,e)=i:^t}gi(z)P t (cos6). (13) 

1=0 * 



As a rule, only the first few harmonics (1=0, 1, . . ., 
L<7) are at one's disposal, but they suffice to specify 
g(z,6). When (12) and (13) are substituted into (4), 
it follows from the orthogonality of the Legendre 
polynomials that 



» i 2Z + 1 2n+l 

1=0 n=Q 



2 



J P l (cosd)P n (cos6)sm6dd (14) 



A2/ + 1 -,. 
1=0 * 



Thus, one can limit oneself to the calculation of the 
first few response function expansion coefficients 



<*i= J iK0)P*(cos0)sin 



1=0, 1, 



L. (15) 



With these coefficients, one is in a position to evaluate 
G(z) for any angular distribution g(zfl) through 
simple superposition according to (14). 

Numerical evaluations of G have been carried out 
for detectors in various positions in box-shaped 
shelters, on the basis of the angular distribution 
functions shown in figure 2. These results apply as 
follows : 



Problem 


Angular 
distribu- 
tion 


Infinite medium 

dose-rate 

rads/hour 

curies/foot 2 


Underground shelter or 
house with roof covered 
by fallout, roof thick- 
ness /i z = 4(0.896 lb of 
concrete per square 
inch) . 

Isotropic comparison dis- 
tribution. 

Sky shine (scattered radia- 
tion) into open hole. 


(a) 

(b) 
(c) 


/(z) =0.983 


jf.00=33.5 



The results are presented as functions of the 
following dimensionless parameters: 



Horizontal eccentricity = 



width _ 
length 



Vertical eccentricity; 



_Vroof-area_V € _ 
height h 



(16) 



(17) 



The range of parameters considered includes < e < 1 
and 0.2<t<20. This appears sufficient for all 
practical applications. 
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Figure 7. Geometric protection factor versus vertical eccen- 
tricity, for fixed horizontal eccentricity e = l. 

The curves labeled (a), (b), and (c) are based on the use of the corresponding 
angular distributions in figure 2. The triangles, solid circles, and hollow circles 
represent results obtained with the response function (7) for a cylindrical 

shelter. 



Figure 7 contains plots of G versus r, for a shelter 
with a square cross section (e=l). Figure 8 con- 
tains plots of G versus e, for different fixed values of 
r. Results are shown for the detector in a corner, 
and in a center position (a=j8=-|). The figures 
confirm quantitatively what one would expect on 
the basis of a commonsense estimate. The pro- 
tection provided by the shelter increases rapidly as 
r is increased, but is generally a very slowly varying 
function of e for fixed r, except when <= has a vftiy 
small value. Dose rates in the corner are on the 
order of four times smaller than dose rates in the 
center. The angular distribution (a) with its strong 
forward peak gives rise to ^-values 2 to 3 times 
larger than those obtained with the diffuse sky- 
shine distribution (c), whereas there is little difference 
between cases (c) and (b). 

The question may be asked to what extent the 
results for the geometrical protection factor are 
affected by the choice of a rectangular shelter cross 
section. To obtain a partial answer, calculations 
were carried out for a detector at the center of a 
cylindrical shelter, with the use of the response 
function (7). As can be seen from these values 
plotted in figure 7, the results thus obtained are, for 
the same value of r, practically identical with the 
results for a detector at the center of a box-shaped 
shelter with 6=1. This means that the geometric 
protection factor depends very little on the shape 
of the cross section. Combined with the previously 
demonstrated insensitivity of G with respect to the 
horizontal eccentricity e, this has the following 
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Figuke 8. Geometric protection factor versus horizontal eccen- 
tricity, for various constant values of the vertical eccentricity. 

The curves labeled (a), (b), and (c) pertain to the corresponding angular dis- 
tributions of figure 2. (A) Detector in center; (B), detector in corner. 



implications: (1) The geometric protection factor 
at the center of a shelter of any shape is, in good 
approximation, only a function of the vertical 
eccentricity r; (2) under these circumstances, a 
calculation using the very simple response function 
(7) will be adequate. 

Center and corner values provide upper and lower 
limits of G. To indicate the transition between the 
two extremes, a systematic mapping of G is presented 
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in table 1 for numerous detector positions in the 
horizontal plane, the shelter being box-shaped with 
e=l and r=2. 

Table 1. Geometric protection factor G as a function of the 
position in the horizontal plane, for a detector in a shelter 
with horizontal eccentricity e = l and vertical eccentricity 
t = 2, calculated for the angular distribution (a) of figure 2 

The position indicators a and are those defined in figure 6. Note that the 
ransformations a->l— a, fi->l—fi, or (a->j8, 0-»a) leave G unchanged. 
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6. Corrections, Accuracy 

6.1. Penetration Through the Side Walls 

This effect has, so far, been investigated only in a 
limited manner to obtain the order of magnitude of 
the necessary correction. To simplify the calcula- 
tions, they were done for a detector at the center of 
the bottom of a cylindrical cavity. For the case of 
an underground shelter (the problem in section 2.1), 
the response function (7) could be amended as folio ws: 






-,<e<T, m=o 



(70 



where 



z'=H—R cot 



The meaning of (7') is that the surface elements of 
the side wall contribute to the dose in proportion to 
the intensity and angular distribution of the radi- 
ation at the distance z—z f of the surface elements 
from the source plane. 

Sample calculations have been carried out for a 
shelter with roof thickness juo2=4, horizontal eccen- 
tricity 6=1, and a depth of 8 ft. (The side-wall cor- 
rection differs from the uncorrected response func- 
tion in that it depends on the absolute size of the 



shelter.) It was found that for vertical eccentricities 

t=1, 2, and 4, 

the increase of the dose is 7, 3, and 1 percent, respec- 
tively. 

6.2. Backscattering 

Some radiation will be scattered toward the detec- 
tor by reflection from the side walls, floor, and ceiling 
of the shelter. The order of magnitude of this effect 
can be estimated roughly on the basis of available 
Monte Carlo calculations of backscattering. The 
results of reference [6] suggest a dose increase by 
approximately 15 percent. 

6.3. Over-all Correction and Error Estimate 

In the final equation for the dose received by the 
partially shielded detector, 



D(z) = Cf(z)G(z), 



(3) 



the factor C, incorporating side-wall and backscat- 
tering effects, is estimated to be 1.2. It is further 
estimated that the factors in (3) have the relative 
errors (SfO)/fG^ ±0.2 and («7)/<7«±0.2, which 
leads to the following error estimate for D: 



§=V(f )*+(£)*- *«• 



(18) 



The authors are indebted to L. V. Spencer for 
valuable discussions, and to J. Hubbell for help with 
the calculations. 
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